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We describe the experimental implementation of a superluminal {i.e. faster than light in vacuo) 
polarization current distribution that both oscillates and undergoes centripetal acceleration. Theo- 
retical treatments predict that the radiation emitted by each volume element of the superluminally- 
moving distribution pattern will comprise a Cerenkov-like envelope with two sheets that meet along 
a cusp. Correspondingly, the emission from the experimental machine is found to be tightly beamed 
in both the azimuthal and polar directions. The beaming is frequency independent and has a 
sharply-defined and unchanging geometry determined only by the speed and path of the moving 
distribution pattern, i.e. by the parameters governing the structure of the Cerenkov-like envelopes. 
In addition, over a restricted range of angles, we detect the presence of cusps in the emitted radi- 
ation. These, which are due to the focusing of wave fronts on a propagating space curve, result in 
the reception, during a short time period, of radiation emitted over a considerably longer period 
of (retarded) source time. The intensity of the radiation at these angles was observed to decline 
more slowly with increasing distance from the source than would the emission from a conventional 
antenna. The angular distribution of the emitted radiation and the properties associated with the 
cusps are in good quantitative agreement with theoretical models of superluminal sources once the 
effect of reflections from the earth's surface are taken into account. In particular, the prediction that 
the beaming and the slow decay should extend into the far zone has been tested to several hundred 
Fresnel distances (Rayleigh ranges). The excellent agreement between the theoretical calculations 
and the data suggests that the apparatus achieves precise and reproducible control of the polariza- 
tion current and that similar machines could be of general interest for studying and utilizing the 
novel effects associated with superluminal electrodynamics. 

I. INTRODUCTION 

Serious scientific interest in the emission of electromagnetic radiation by charged particles moving faster than the 
speed of light in vacuo began with the work of Sommerfeld in 1904 and 1905 The timing of this work was perhaps 
unfortunate, as the special theory of relativity was published within a few months |2(; any known particle that has 
a charge has also a rest mass and so cannot be accelerated through the speed of light. Moreover, an additional 
conceptual barrier discouraged further work: even a massless particle cannot move faster than light in vacuo if it is 
charged, because if it did, it would give rise to an infinitely strong electromagnetic field on the envelope of the wave 
fronts that emanate from it 

By and large, this field remained dormant until Ginzburg and colleagues 0, 0, pointed out that, though no 
superluminal source of electromagnetic fields can be point-like, there are no physical principles preventing extended 
faster-than-light sources. The coordinated motion of aggregates of subluminally-moving charged particles can give 
rise to macroscopic polarization currents whose distribution patterns move superluminally 0, El 0] ■ Further relevant 
progress occurred with the theoretical prediction that extended sources that move faster than their own waves could be 
responsible for the extreme properties of both the electromagnetic emission from pulsars (rapidly spinning, magnetized 
neutron stars) 0, 0, H, El El an d the acoustic emission by supersonic rotors and propellers ■ Finally, theoretical 
treatments of a particular class of source of this type [l2L IT3 | were stimulated by the design and construction of the 
apparatus that is the subject of the current paper. This apparatus implements an extended polarization current 
distribution that has the time dependence of a traveling wave with centripetally-accelerated superluminal motion and 
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an oscillating amplitude. Refs. point out several distinctive characteristics of the emission by such sources: 

1. The radiation emitted from each volume element of the superluminally- moving polarization-current distribution 
pattern comprises a Cerenkov-like envelope with two sheets that meet along a cusp. The combined effect of the 
volume elements leads to a frequency-independent beaming in both the azimuthal and polar directions whose 
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geometry is determined by the structure of the Cerenkov-like envelopes 0, IT2I Il3l ] . 

2. The above-mentioned cusps, which spiral upwards and outwards from the source into the far zone |6L Il2l Il3| . 
result from the centripetal acceleration of the source; without acceleration, the wave fronts would merely pile up 
on a two-dimensional envelope [|| |nj . The cusps result in the reception, during a short time period, of radiation 
emitted over a considerably longer period of (retarded) source time. Consequently, the intensity of the radiation 
at these angles declines more slowly with increasing distance from the source than would the emission from a 
conventional antenna. Note, however, that energy conservation is not violated; large intensities on the cusp are 
compensated by weaker radiation fields elsewhere. The propagating cusp is being reconstructed constantly from 
conventional (i.e. spherically-decaying) waves that combine and then disperse. 

3. Under certain conditions, the emitted radiation contains a broad spectral distribution of high frequencies that 
are not present in the synthesis of the moving, oscillating polarization current. These originate from the 
centripetal acceleration, and from space-time discontinuities forced on the polarization-current distribution by 
the constraints of its circular path. 

Based on the above ideas, we have constructed an experimental realization of an extended, oscillating superluminal 
source that undergoes centripetal acceleration. The present paper describes the basic principles of the machine, and 
its use to demonstrate the first two theoretical findings listed above; the high frequency emission will be the subject of 
a future work. We find that the emission of the source is well described by the models of Refs. 0,E1; confirming the 
validity of extended superluminal polarization currents as precisely controllable sources of electromagnetic radiation 
with novel properties. 

The paper is organized as follows. Section ^ gives a description of the experimental considerations necessary to 
understand the data from the superluminal source. Since practical superluminal sources are relatively novel, and 
since the raw experimental data are affected by the proximity of the source to the ground in a way familar only to 
communications engineers ^ij , this section is by necessity rather detailed; without a thorough explanation of these 
effects, the experimental data are not readily understood. Within it, Subsection III Al describes how the polarization 
current distribution is created and animated (set in motion) at superluminal speeds within a curved strip of alumina. 
The experimental geometry, and in particular its relationship to the coordinate systems used in the theoretical 
papers 0, ^| is discussed in Subsection III Bl whilst the detection system and the effects of interference from the 
ground are described in Subsection III CI Experimental results commence in Section IIIII which describes data that 
conclusively demonstrate the three-dimensional angular distributions of the Cerenkov-like envelopes emitted by the 
source. Section |W] describes range tests exhibiting the non-spherically-decaying radiation close to the cusp; in spite of 
the non-ideal nature of the experimental machine and the measurement geometry (which complicated by interference 
from the ground) it is possible to conclusively identify this feature within the emission data. A discussion, conclusions 
and summary are given in Section E] The paper ends with an Appendix describing the mathematical treatment of 
interference due to reflections from the ground. In addition, the Appendix details the fitting methods used to obtain 
the exponent of the non-spherical decay associated with the cusps, and summarizes the parts of the theory given in 
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Refs. 0, that are used to model the beaming observed in the experimental data. 

II. EXPERIMENTAL DETAILS 

A. Animation of a superluminal polarization current 

The first step in investigating superluminal effects is to identify a source that can travel faster than the speed of 
light in vacuo without violating Einstein's theory of special relativity . The required source explicitly appears in 
the Ampere- Maxwell equation, which underlies the emission of electromagnetic radiation |l5lll6j |: 

^ , __ _ 9D <9E dP ,_ T . , 

V x H = J frcc + — = J free + e Q — + — (SI units). (1) 
at at at 

Here H is the magnetic field strength, Jf ree is the current density of free charges, D is the displacement, E is the 
electric field and P is the polarization; the term dP/dt represents the polarization current density. An oscillating or 
accelerating current density Jf ree emits electromagnetic radiation (i.e. it acts as a source) and is well-understood; this 
is the basis of conventional radio transmission 0,01 an d the use of synchrotrons as light sources 0|. However, 
as a current density of free charges is composed of particles such as electrons or protons, Einstein's theory of special 
relativity naturally forbids such a source from traveling faster than the speed of light. By contrast, there is no such 
restriction on a polarization current. Fig. Ea)-(b) shows a simplified dielectric solid containing negative and positive 
ions. In (b), a spatially- varying electric field has been applied, causing the positive and negative ions to move in 
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FIG. 1: Left: experimental animation of a superluminal polarization current, (a) A simplified dielectric solid containing negative 
(©) and positive (©) ions. In (b), a spatially- varying electric field has been applied, causing the positive and negative ions to 
move in opposite directions; a finite polarization P has therefore been induced. If the spatially-varying field is made to move 
along the direction of the arrow, the polarized region moves with it. (c) Schematic side view of a practical superluminal source, 
showing metal electrodes above a strip of dielectric (shaded region) and a ground plate below it. "0" indicates that there is 
no voltage on that particular upper electrode; the symbol + indicates a positive voltage applied to the upper electrode. The 
voltage on the electrodes produces a finite polarization of the dielectric (darker shading) . (d) By switching the voltages on the 
electrodes on and off, the polarized region (darker shading) can be made to move along the dielectric, (e) Top view, showing the 
curvature of the dielectric (lighter shaded region). The curvature introduces centripetal acceleration in the moving polarized 
region. Note also that the electrodes (black shading) cover only part of the top surface of the dielectric. In the discussion that 
follows, we refer to this section of the source (electrodes, dielectric and groundplate) as the "array", for convenience. Right: a 
photograph of the experimental array; the copper shields have been removed from the 41 amplifiers so that their components 
are visible. 



Symbol 


Definition 


u> 


Angular rotation frequency of the distribution pattern of the source. 


n 


The angular frequency with which the source oscillates (in addition to moving). 
One of the two angular frequencies used to drive the experimental source (Eq.^J. 


f, f± 


Frequency of the radiation generated by the source; f± = \ fl ± mui\/2n. 


n — 2nf /lu 


The harmonic number associated with a particular radiation frequency. 


m 


The number of cycles of the sinusoidal wave train representing the azimuthal 
dependence of the rotating source distribution [see Eq|5| around the 
circumference of a circle centred on, and normal to, the rotation axis. 


r\ = tjiu 


The second of the two angular frequencies used to drive the experimental source (Eq[2j. 



TABLE I: Definition of the various frequencies and numbers used to describe the source and the emitted radiation, after 
Refs. [H[ll. 

opposite directions. A finite polarization P has therefore been induced. If the spatially-varying field is made to move, 
the polarized region moves with it; we have a traveling "wave" of P (and also, by virtue of the time dependence 
imposed by movement, a traveling wave of dP/dt). Note that this "wave" can move arbitrarily fast (i.e. faster than 
the speed of light) because the individual ions suffer only small displacements perpendicular to the direction of the 
wave and therefore do not themselves move faster than the speed of light. 

Thus far, the bulk of the theoretical analysis of superluminal light sources has been carried out for rotating, oscil- 
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lating polarization current distributions |l2l , which possess a polarization given (in cylindrical polar coordinates) 
by 

P(r, <p, z, t) = s(r, z) cos[m((/? — wi)] cos tit. (2) 

Here s is a vector field describing the orientation of the polarization; it vanishes outside the active volume of the 
source |l2j . The other parameters are defined in Table [I] The experimental superluminal source described in this 
paper represents a discretized version of Eq. [21 It consists of a continuous strip of alumina (i.e. the material to 
be polarized by the applied electric field) on top of which is placed an array of metal electrodes; underneath is a 
continuous ground plate. This forms what is in effect a series of capacitors (a schematic is shown in Fig. Uc)-(e)). 
Each upper "capacitor" electrode is connected to an individual amplifier; turning the amplifiers on and off in sequence 
moves a polarized region along the dielectric (Fig.^c)-(d)). 

In practice, the dielectric in the experimental source is a strip corresponding to a 10° arc of a circle of average 
radius a = 10.025 m, made from alumina 10 mm thick and 50 mm across. Above the alumina strip, there are 41 
upper electrodes of mean width 42.6 mm, with centers 44.6 mm apart. Whilst the ground plate is the same width as 
the dielectric strip, each electrode covers only the inner 10 mm of the upper surface (shown schematically in Fig^e)). 
This has the effect of producing a polarization in the alumina with both radial and vertical components (see Eq. EJ . 
In the discussion that follows, we shall refer to the assembly consisting of groundplate, dielectric and electrodes as 
"the array" for convenience. 

Each upper electrode of the array is driven by an individual shielded amplifier just behind it; to simulate Eq.[2 the 
jth (j — 1, 2, 3....) electrode is supplied with a voltage 

Vj = Vq cos[r]{jAt - t)] cos fit. (3) 

Comparison with Eq. [2] shows that the angular frequency r\ = mui; At = Atp/ui, where Aip — n/(18 x 40) is the 
angle subtended by the effective center separation of adjacent electrodes. The velocity v with which the polarization 
current distribution propagates along the dielectric is set by adjusting At to give v = aAif/At (see Fig. E{a)). Given 
the dimensions of the experimental array, v > c is achieved for At < 148.8 ps. (At this point it should be stated that 
the modulation scheme used to animate the superluminal polarization current is quite distinct from any used thus far 
in conventional phased-array antennae |19|.') 

The application of the electrode voltages produces a stepped approximation to a sinusoidal polarization-current wave 
that moves at speed v along the dielectric (see Fig. Ufa)). Detailed analysis of this type of waveform (see Ref. [l^. 
Appendix) shows that the dominant emission is almost identical to that produced by a pure sinusoidal wave with the 
same wavelength and amplitude. Most of the emission occurs at two frequencies, f± = | d= 77 1 /2tt = |fi ± mcj|/27T, as 
indicated in Table [I] 

As mentioned above, the speed of the source (i.e. the speed at which the sinusoidal polarization current "wave" 
propagates along the dielectric) is set by the parameter At. Experimentally, this is done by observing the time 
dependence of the voltages on each electrode (Figure Efb)) using a high-speed oscilloscope. 

Finally, it should be emphasised that great care was taken to minimize unwanted radiation from the amplifiers, 
cabling and signal generators of the apparatus to ensure that the detected signal was dominated by the emission 
from the polarization current. Steps taken included proper encasing of all of the signal generation equipment, the 
enclosure of the 41 amplifiers inside thick copper shields (not shown in Fig. ^| and the correct impedance matching 
of connections. 



B. The experimental geometry 

To map out the spatial distribution of the radiation from the superluminal source, one would ideally position the 
array in empty space and move a detector to points around it, covering a comprehensive range of angles and distances. 
However, practical experiments must be conducted close to the earth's surface, where reflections from the ground are 
a complicating factor [14| (see Appendix); numerous tests showed that highly reproducible data are obtained when 
the "ground" under the path between array and detector is smooth, level and of a uniform quality 20] . After trials at 
several locations, the bituminised macadam surface of an airfield runway (Turweston Aerodrome, near Northampton, 
UK) was found to be almost ideal in providing a substantial length of well-characterized surface. All of the experiments 
reported in the current paper employ the runway, which is approximately 900 m long and 20 m wide. A change in 
slope about two-thirds of the way along the runway restricts the range tests reported in Section IIVI to distances of 
600 m or less. However, measurements of the angular distribution of the radiated power from the array (Section ^Hj 
employed the full length of 900 m. 
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FIG. 2: (a) The upper curve shows the instantaneous values of the voltages on the 41 upper electrodes of the array; r]/2Ti = 
552.653 MHz and At = 140 ps, corresponding to v « 1.064c for a = 10.025 m and Aip = tt/(18 x 40) (see Eq. The 
resulting polarizations under the electrodes will be proportional to these voltages. The lower curve (offset for clarity) gives the 
instantaneous voltages 420 ps later; the polarization "wave" has moved along the dielectric. (Note that the cos fit term in Eq.0 
hardly changes during this time.) (b) Example of the time-dependent voltages on array electrodes 20 and 21 (upper curves) 
and 20 and 26 (lower curves) for the same parameters as in (a). Note that the lower- frequency temporal modulation of each 
electrode voltage produced by the second cosine term in Eq. [3] (dotted curves) is in phase for all electrodes. 



The geometry of the runway dictates that the best configuration for measuring the spatial and angular distribution 
of the radiation is to fix the position of the detector (on the runway, at a distance R from the source) and to change 
the orientation of the array to access a full range of polar and azimuthal angles. This is accomplished by placing the 
array on a platform that tilts about a horizontal-axis pivot at one end; this is in turn fixed on a turntable that can 
rotate about a vertical axis. 

It is possible to mount the array on the platform in two configurations. In the first, the array plane is vertical 
(Fig, ffia)). so that its long axis can be tilted (about the horizontal pivot) to an angle eV with respect to horizontal; 
the turntable is then rotated about a vertical axis so that the array plane makes an angle t9y (i n the horizontal 
plane) with the line connecting the detector at P and the array center. In the second configuration, the array plane 
is initially horizontal (Fig. |3Jb)); the array plane is then tilted using the horizontal pivot to an angle $h with respect 
to horizontal. The turntable is then rotated about a vertical axis so that the projection of the long axis of the array 
on the horizontal plane makes an angle (f>n with the line connecting the detector at P and the array center. 

Finally, to aid comparison of experiment and theory, it is necessary to relate the experimental spherical polar 
coordinates (R, iVj^v) and (-R, $ti, 4>n), which are measured relative to the array center, to the cylindrical and 
spherical polar coordinates used to evaluate the integrals in the Appendix and in Refs. |12,ll3j. In the latter papers, 
each volume element of the superluminal source distribution traverses a complete circular orbit centered at the origin 
of the coordinate system, and the radial distances to the individual source elements and to the observation point at 
(i?p, 9p, (fip) are specified with respect to the common center of these circular orbits (see Fig.^c) and (d)). Although 
the experimental source does not encompass a complete circle, the high degree of symmetry represented by this 
coordinate system greatly facilitates the computation of the integrals in the Appendix. The relationship between the 
experimental coordinates and the coordinate system used in the Appendix is shown in Figs.|3Jc) and (d). 
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FIG. 3: Left: the coordinates used to measure the spatial distribution of radiation from the array. In the experiments ((a),(b)), 
the detector is at a fixed point P, distance R from the array; the array is placed on a platform tilting about a horizontal axis 
pivot that is in turn fixed to a turntable rotating about a vertical axis. In (a) , the array plane is vertical; its long axis is tilted 
to an angle 0v with respect to horizontal. The turntable is rotated about a vertical axis so that the array plane makes an angle 
#v (in the horizontal plane) with the line connecting P and the array center. In (b), the array plane is initially horizontal; 
the array plane is then tilted to an angle i?h with respect to horizontal. The turntable is rotated about the vertical axis so 
that the projection of the long axis of the array on the horizontal plane makes an angle (j>u with the line connecting P and 
the array center. Note that the coordinate systems are not equivalent; tan^v = tan </>h/ cos i?h and sini?v = sin i?h cos 0h • 
(c) The array and detector positions in terms of the coordinate system used in the Appendix and Refs. [T^. Il3|. In the model, 
each volume element of the sourse moves on a circular path centered at the origin O; a typical path is shown, with the array 
being represented by thicker shading. The array center is at (a, 0,0), and the tangent to the array is parallel to the y axis; 
the distances, from O, to the detector (Rp) and to the array center (R) are compared, (d) Definition of the spherical polar 
coordinates Rp,9p,ipp used in the Appendix. The arrow indicates the direction of propagation of the polarization current 
distributions in the experimental array. Right: the array in use during the experiments. 



C. Detector system, array and detector heights, and interference from ground reflections 

The detector was a Schwarzbeck Mess-Elektronik UHA 9105 half-wave tuned dipole antenna, optimized for the 
frequency range 300 MHz to 1 GHz. and mounted on a mobile wooden tower that allowed it to be fixed at varying 
heights hr> above the ground. The antenna was fitted with an internal coaxial BAL-UN to ensure good impedance 
matching with an external load and was connected via a 6 m long impedance-matched coaxial cable (50 f2) to an 
Anritsu MS2711B portable spectrum analyser. The background noise floor with this equipment was between -120 
and -130 dBm, depending on weather conditions at the airfield. 

In the current paper, the settings of the experimental machine are such that the emitted radiation is polarized 
chiefly in directions parallel to the plane of the array. When the array is orientated as in Fig. [3Ja) , the detector aerial 
(see below) is arranged to detect vertically-polarized radiation; if the array is orientated as in Fig. E^b), the detector 
is configured to detect horizontally-polarized radiation. 

For the measurements described in the following sections, the array turntable was mounted on a scissor lift that 
allowed the height of the center of the array above the runway to be varied between /ia = 3.3 m and 5.0 m; for each 
measurement, the heights /ia and hr> were measured to an accuracy of ±5 cm. An accurate measurement of /ia and 
/id is important because the intensity received at the detector results from interference between signals traveling along 
a direct path between array and detector and reflections from the ground [l4|. The situation is shown in Fig. Ufa). 
This process is covered in detail in the Appendix. 

The effect of the interference on the signal received from a simple horizontally-polarized dipole aerial is shown in 
Fig Eft). Firstly, clear fringes are seen in the detected intensity as a function of distance. A second consequence of 
the interference is less obvious. Fig. 01 shows received power in the logarithmic units Belm versus log 10 (distance); the 
gradient of the plot should be —7, where (received power) tx (distance)" 7 . Naive expectations based on the emission 
of a dipole surrounded by infinite empty space lead one to predict 7 = 2 (the inverse square law [T^l: in practice, 
the interference leads to 7 rj 4 (see Fig. 0Jb)). Such a dependence is familiar to radio engineers as the penultimate 
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FIG. 4: (a) Geometry of the interference between a direct wave from the source S, height /ia above the ground, and a wave 
reflected from the ground. The detector, at height /id above the ground, is labeled by P. The direct waves travel a path length 
denoted by R. The path source-ground-detector is equal in length to a direct path from the source to an image of P at P'. The 
length of this equivalent path is denoted by R T . The other symbols shown are used in the Appendix, (b) Effect of interference 
on the power received from a horizontally-polarized dipole aerial excited at 598.696 MHz as a function of distance (h&. = 5.12 m; 
/id = 5.0 m). The points are data, and the straight line is a least-squares fit to the data beyond the final interference fringe; 
the gradient is 7 = —3.8. 

term of the "Egli" relationship |2lJ describing the power lost for frequencies / in the hundreds of MHz range when 
transmitting over a path close to the earth's surface: 

path loss (dB) = 113 + 201og 10 / + 401og 10 (path length (miles)) - 20 log 10 (/i A /iD). (4) 

The effects of interference obviously complicate the search for the non-spherical decay associated with the cusp (see 
Section llVjl . One criterion used in the experiments was to look for significantly smaller path losses than Egli (Eq.^J; 
more detailed analysis was then used to model the detailed behavior of the data. 

Reflections from the ground complicate also the modeling of the angular distribution of the emission from the array. 
A detailed treatment of the coordinates needed for delineating the direct and reflected rays is given in the Appendix, 
as are the equations from Refs. mm required for modeling the data. 

III. EXPERIMENTAL RESULTS: INTENSITY DISTRIBUTION OF THE RADIATION 

A. Intensity as a function of $v at fixed 4>v 

The array was mounted in the configuration shown in Fig. Eta) and the detector antenna placed at a fixed distance 
R away; detector and array were fixed at heights /id and /ia respectively above the level runway. The angle (f>y was 
set to a series of fixed values with an accuracy of ±0.5° using a Suunto clinometer, and then the angle $v was varied 
using the turntable. The measurement of ??v w & s made to an accuracy of ±0.5° using a graticule on the turntable. At 
each value of (i?V; <^v); the intensity of the signal at the detector antenna was measured using the spectrum analyser. 

Fig.[Sfa) shows the effect of varying the source speed v = aAip/At on the intensity distribution as a function of #v- 
Data for v/c = 1.064, 1.25 and 2.0 are shown; for each data set, ?y/27r = 552.654 MHz, Q/2tt = 46.042 MHz, and the 
observation is made at the higher of the two emitted frequencies, /+ = |(0 + mo;)|/27r = (fi + r])/2ir = 598.696 MHz. 
As the source speed increases, the peak in detected power moves out to higher values of fly, in line with expectations; 
the contributions that are received from different source elements, by an observer in the far field, are in phase when 
(see Appendix, Eas. 1251 and l4"7|l 

*}• (5) 



arcsm I . 

V 77,7; / 



or t?v = arcsin 



1 - 



/mcy 
V 77,7; / 



Here n = 2nf/uj is the harmonic number of the emitted radiation with / = / + or /_, and m — rj/uj is determined by 
the angular frequency u> — v/a — Aip/At of the rotation of the source (see Table [IJ 12]. Note that the data are well 
reproduced by the model (see Appendix, Eq. I44fl using the experimental source speeds. 
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FIG. 5: (a) The effect of varying the source speed, v = aAip/ At on the angular distribution of the radiation from the array. 
For each data set, <f>\ — 2° and the array to detector distance is R — 600 m. Data are points, identified as follows: filled circles: 
v — 1.064c, /ia = 3.7 m, = 4.1 m; crosses: v = 1.25c, Ha = ho = 5.0 m; diamonds: v — 2.0c, Ha = ho = 5.0 m. The curves 
are fits, using the model (Eq. described in the Appendix. For each data set, rj/2-K = 552.654 MHz, Q/2tt = 46.042 MHz, 
and the observation is made at the higher of the two emitted frequencies, /+ = |(fl + mu>)\/2ir = (f2 + n) /2n — 598.696 MHz. 
Note that the peak emission moves out to higher values of $v as the source speed increases, (b) Power detected versus f?v over 
a wider angular range than (a); here, R = 400 m, /ia = ho = 5.0 m, <f>y = 2° and the frequency is the same as in (a). The 
points (hollow diamonds) are data; the model (using the experimental source speed) is the continuous curve. 



Compared to the emission of a conventional antenna |25j , the rather small side lobes are a further characteristic of the 
beaming produced by an accelerated superluminal source. To show this more clearly, a single data set, corresponding 
to v/c — 2.0, is shown in Fig. EJV). Although the projection of the length of the array onto a plane normal to the 
observation direction is comparable to or smaller than the radiation wavelength and the observation is made at many 
hundreds of Fresnel distances, most of the power is transmitted in a beam around 20° wide, with small side lobes. 
This beam structure is reproduced very well by the model (Appendix, Eq. I44fl , using the experimental source speed as 
an input parameter. (In the case of the present emission mechanism, the side lobes are caused mainly by interference 
from the ground; it follows from Eq. [33] that in free space they would be much weaker.) 

Whilst Fig. [5] shows data for the emission at /+ = |(f2 + r])\/2n, the description of the source, Eq. |5] leads one 
to expect simultaneous emission at a second frequency, /_ = |(fi — r})\/2ir. As /+ and /_ have different harmonic 
numbers, n, their angle of emission will differ slightly (see Eq. |S]and Table P). This effect is shown for source speeds 
of v/c = 2.0 in Fig. HJa) and v/c — 1.25 in Fig.HJb). In each case, the emission at /+ is peaked at a higher value of 
$v than that at /_, as predicted by Eq. |S]and the solid curves (see Appendix). 

It has already been mentioned in the discussion of Fig. [SJa) that the beam is rather sharply defined as a function 
of $Vj even though the projection of the length of the array on the plane perpendicular to the observation direction 
is comparable to or smaller than the radiation wavelength. This characteristic is independent of distance over the 
range measured, as shown in Fig. [7| which plots data (points) for 200 m < R < 900 m. For all of these R values, 
each representing hundreds of Fresnel distances, the central part of the beam preserves its overall width and general 
shape. This development of beam structure with R is reproduced by the theoretical description (Fig. [7| curves; see 
Appendix, Eq. I44|) , using the experimental source speed as input parameter. 

Note that all of the data sets plotted in this Section involve small values of 0v(~ 2°). The reason for this constraint 
is that as well as having a peaked structure as a function of the polar angle t?Vj the emission of the source is also 
peaked in the plane of rotation (i.e. around 0y = 4>n — 0), as expected in the theoretical description 0,0. This 
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FIG. 6: Comparison of the emission at /+ = (77 + n)/2yr = 598.696 MHz (hollow diamonds) or /_ = (j]-Q)/2tt = 506.614 MHz 
(filled circles) for two different source speeds. In (a), v = 2.0c, R — 400 m, h\ ~ 3.32 m, ho = 3.71 m and 0v = 2.5°. In (b), 
v = 1.25c, h\ = ho = 5.0 m, <j>v = 2° and R — 600 m. For all data sets, the values of Q and rj are the same as those in Fig. 
and in all cases the curves are generated by the model fEa. 144^ using the experimental values of f± and v/c. Note that, in 
agreement with Eq. the emission is peaked at higher values of i?v for /+ than for /_ . 
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FIG. 7: (a) The #v dependence of the detected power for R — 200 m (circles), 500 m (triangles), and 900 m (diamonds), 
(b) The $v dependence of the detected power for R — 300 m (squares) and 600 m (inverted triangles). In both cases, 
/+ = (rj + n)/27r = 598.696 MHz, v = 2.0c, fiA = ho = 5.0 m and 0v = 2°; points are data and curves are fits to the model in 
the Appendix (Eq. . 
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FIG. 8: The (f>n dependence of the detected power (logarithmic scale) for R — 200 m (upper) and 900 m (lower) for #h = 30° 
(filled circles) , 35° (hollow diamonds) and 40° (inverted triangles) . Solid curves are fits to the model in the Appendix (Eq. 1441 . 
The frequency is /+ = (77 + Q.)/2ty = 598.696 MHz, v = 1.064c, h A = hu = 5.0 m. 



observation will be described in the following Section. 



B. Intensity as a function of 4>h at fixed $h 

The array was mounted in the configuration shown in Fig. GJb) and the detector was placed at a fixed distance R 
away; detector and array were fixed at the respective heights of /id and /ia above the level runway. The angle $h was 
set to a series of fixed values to an accuracy of ±0.5° using a Suunto clinometer, and then the angle <fin was varied 
using the turntable. The measurement of 4>~n was carried out using a graticule on the turntable to an accuracy of 
±0.5°. At each value of (tin, <^h) the intensity was recorded. 

Fig. IS] shows the beam profile as a function of 0h for $h = 30°, 35° and 40° and distances R = 200 m and 900 m. 
Note that the peak emission is centered around </>h = 0, as expected from theory [l^; variations of </>h of around ±10° 
away from zero do not affect the intensity much, but thereafter, the intensity falls away rapidly. 

Taken in conjunction with Fig. the data in Fig. [H] show that the emission from the super luminal source is tightly 
beamed in both the azimuthal and polar directions. The observation of a peaked emitted power distribution (see 
Fig. HJJ) at two frequencies shows that this effect is frequency independent, depending instead on the speed of the 
source and the geometry of the array. Once again we contrast this with the case of a conventional antenna j2f| ; the 
beaming occurs even though the projection of the length of the array onto a plane normal to the observation direction 
is comparable to or smaller than the radiation wavelength and the observation is made at many hundreds of Fresnel 
distances. 

Finally, we remark that the evolution of the beam profile as $h and vary is in good agreement with the 
expectations of theory 0| ■ Fig- E9 shows the the 0h dependence of the detected power (linear scale) for several 
values of z?h- The evolution of the beamshape, including the development of a precursor (a side maximum), is well 
reproduced by the model (Appendix, Eg . I44fl . using the experimental value of v/c. 
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FIG. 9: (a) The <j>n dependence of the detected power (linear scale) for i?h = 40° (hollow circles), 35° (filled diamonds), 30° 
(filled squares) and 25° (filled circles). Points are data; the curves are the model (Eg. 1441 predictions, (b) The 0h dependence 
of the detected power for i?h = 20° (inverted triangles), 15° (diamonds) and 10° (crosses). Points are data; the curves are the 
model predictions (Eq.|HJ. In both (a) and (b) R = 400 m, the frequency is /+ = (?? + Q)/2tx = 598.696 MHz, v = 1.06c, 
h,A = ho = 5.0 m. 

IV. NONSPHERICALLY DECAYING COMPONENT OF THE RADIATION 

A. General properties 

The previous section has shown that the angular distribution of the emission from the experimental source is in good 
quantitative agreement with theoretical predictions for a rotating, oscilla ting superluminal source (see Appendix) [T^ . 
Il3j . However, the most controversial aspect of the model of Refs. |ol IT^. Il3| has been the prediction of a component 
of the emission that decays nonspherically, i.e. as 1/i?, rather than 1/R 2 This component is due to an effect, 

only possible in the case of an accelerated superluminal source, whereby waves that were emitted over a finite interval 
of emission time arrive at the detector during a much shorter interval of reception time. The constructive interference 
of such superposed waves is strongest if their emission time interval is centered at the instant at which their source 
approaches the detector (along the radiation direction) with the wave speed (the speed of light) and zero acceleration. 
The locus of the series of points for which this occurs for a single volume element of the source constitutes the 
propagating space curve which we have here called the cusp 0, Il3| . 

In this section, we describe the location and the characteristic features of this emission for a source speed of 
V = 1.064c and the frequency /+ = (fi + T])/2tt = 598.696 MHz; the theoretical predictions for this speed lead one to 
expect that the cusp loci will be found close to $h = 20° at small values of 0h: For a source that is centered at (p = 
and an observation point that lies in the far field, Eq. 49 of Ref. reduces to <p-p ~ —3tt/2, Op ~ arcsin[c/(raj)], in 
which no = v. 

As a reference point, Fig. HOf a) shows the detected power versus array-to-detector distance R, both plotted on 
logarithmic scales, for the case of the source moving subluminally; the velocity is v/c — 0.875. Note that the data 
(points) are very similar to those shown in Fig. 0] for the emission of a dipole aerial; indeed, the data are well fitted 
using the interference model for a point-like source and spherical decay described in the Appendix (curve), especially 
at larger distances where the finite array size is much smaller than R. This suggests that the emission of the array 
running subluminally is entirely "conventional" ; the detected power obeys the inverse square law (once the effects of 
interference are taken into account), and at large distances the source appears point-like. 

Figure UTir b) shows detected power versus log 10 (i?) for the source running superluminally (v/c — 1.064), for array- 
to-detector distances R measured along the expected direction of the cusp loci. The contrast with Fig. ITUT a) is 
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FIG. 10: (a) Detected power (logarithmic scale) versus log 10 (-R), where R is the array-to-detector distance for a subluminal 
source, velocity v/c — 0.875; Ha = hr> = 5.0 m, </>h = 0, i?h = 35°, /+ = (77 + Q)/2tt = 598.696 MHz. The points are data; 
the curve is the interference model described in the Appendix; the refractive index is N = 1.074. (b) Detected power versus 
log 10 (i?) for the superluminal source close to the expected cusp direction. The source velocity is v/c = 1.064, Ha = hr> = 5.0 m, 
H = -5°, H = 20°, /+ = (»7 + Q)/2tt = 598.696 MHz. Data are points; the dotted curve is a fit to Eq. [TJJ (Appendix). This 
model is valid for larger R and assumes both no beaming and an intensity proportional to 1/R 2k . Fit parameters are 2k = 0.90 
and N — 1.074. The solid curve represents a fit to Eqs. I18l and ll9l that uses a Gaussian beam of width 2w; fit parameters are 
2k = 1.05, 2w — 4.18° and K = 16.6. (c) Detected power versus log 10 (i?) to one side of the cusp direction; all experimental 
details are the same as in (b) except <^h = 0. Data are points, the dotted line is a fit to Eq.[HI(2K = 1.41, N = 1.074). The 
solid curve is a fit to Eqs El and El (2k = 1.47, 2w = 6.69°, K = 17.3). (d) Same as (c) but with ct> H = -10°. Data are 
points; the fit to Eq.[T7|is dotted (2k = 1.47, N = 1.074). the solid curve is a fit to Eqs El and El (2k = 1.48, 2w = 6.56°, 
K — 14.6). (The data shown here are limited to R < 600 m because the runway at the Turweston Aerodrome dips sharply in 
its last 300 m.) 

immediately visible, and involves three distinct differences: 

1. The interference fringe visible at around log 10 (i?) = 2 in (a) is absent in (b). 

2. At larger distances, the detected power falls off with distance more slowly in (b) than in (a). 

3. The data in (b) possess a much greater "scatter" than those in (a). This scatter occurs in spite of the fact that 
the detected power is ~ 50 dB above the noise floor; moreover, the scatter is to a large extent reproducible. 

The last point can be understood by noting that the cusps represent emission by those volume elements of the source 
that travel towards the detector at zero acceleration and the speed of light, volume elements that collectively function 
as a coherent source 01 • The power emitted by a coherent source (such as a laser) normally contains reproducible 
variations known as "speckle" j24[; the reproducible scatter seen in Fig. irUT b) is due to just such a phenomenon. 

Turning to the other differences between Fig. HUf a) and (b), the absence of an interference fringe is expected 
if the radiation along this direction is tightly beamed (the loci of the cusps are expected to have narrow angular 
width 0,EH) j so that it is not reflected from the ground close to the array; a detailed treatment of this effect is given 
in the Appendix (Eas. ITH1 and IT§|) . The slower decrease in power with increasing distance is another consequence of 
the properties of the cusps; we shall return to the exact exponent describing the decrease below. 

As mentioned above, the cusps occupy only a small angular range; morover, they are expected to become narrower, 
the further one goes from the source. Some idea of this effect is gained by rotating the array by 5°; whereas Fig. llOf b) 
has </>h = —5°, Figs. liur c) and (d) show similar data for 0h = and 0h = —10° respectively. Note that the small shift 
in direction has resulted in data that exhibit a smaller degree of "scatter" compared to those in Fig. HUf b) (especially 
at larger i?), and that the detected power falls off more rapidly with increasing distance in (c) and (d) than in (b). 
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FIG. 11: (a)and(b): Ratio P cusp /P sub (see Eq.EJ of the detected power with the source running super luminally (at v/c = 1.064) 
to that with the source running subluminally (v/c = 0.875) versus R, measured on two separate occasions. Both data sets follow 
paths close to the expected direction of the cusps; h\ = ho = 5.0 m, 0h = —5°, $h = 20°, /+ = (77 + fl)/2n — 598.696 MHz. 
Data are points; the lines are fits to Eq.HJwith fj, » 1. A linear fit through the origin implies that the power in the superluminal 
case would decline in free space as 1/R (see Eq. |SJ. (c) and (d) show similar plots just off the cusp direction; (fin = in (c), 
and ifin = —10° in (d). All other parameters are the same as in (a) and (b). The fits (solid curves) are to Eq. HJwith /i = 0.62 
in (c) and /1 = 0.52 in (d); i.e. the power ratio increases more slowly with distance, and cannot be fitted with a line through 
the origin. 



The first observation is entirely expected; once one is off the cusp loci, the source no longer appears as coherent [13J, 
and so the "speckle" effects should decrease. Along directions that appreciably differ from that of the cusp loci, the 
emitted power has the same dependence on distance as in Fig. HUf a). i.e. as in the case of a subluminal source or a 
conventional dipole. 

Having identified the distinct features of the non-spherically-decaying component of the radiation, we describe in 
the following sections three experimental methods for determining the exponent 2k that describes the rate at which 
the intensity of this component decays with distance R (i.e. intensity oc 1/R 2k ). In each case, we shall find that 
2k sa 1, in agreement with theoretical expectations |3. IT^. 



B. Determination of the rate at which power falls off with distance using models of the interference 

Comparison of Figs. HUf a)-(d) shows that the detected power decreases most slowly with increasing distance in 
the case of superluminal velocities and measurements along the cusp direction (Fig. HUf b)). However, obtaining 
the exponent 2k for this decrease is complicated by the effects of interference from the ground; we have already 
seen that even the emission of a simple dipole does not decrease as 1/R 2 (Fig. Efb)) in proximity to the ground. 
The part of the radiation corresponding to the cusps provides a further complication in that it is expected to be 
tightly beamed |fj, [13J ; the current experimental machine agrees with these predictions in that the beam width of 
of the non-spherically-decaying component is appreciably (by an order of magnitude or more) narrower than the 
beams observed in the spherically-decaying part of the radiation (as shown in Figs. I5T§|) . A narrow beam will not be 
appreciably reflected from the ground when the detector is close to the array; hence, features due to interference will 
be suppressed. 

Two treatments for this problem are included in the Appendix (Eas. fT71to ll9l) . If the beam has an angular spread 
2w, then at distances R > hj^/w, a substantial proportion of the beam will be reflected from the ground; therefore 
a relatively simple treatment of interference (Eq. I17|) can be used. The only input parameters are the power law 
of the decay of the radiation with distance (intensity oc 1/R 2k ) and N, the refractive index of the ground (already 
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constrained to be N = 1.074 by fitting the interference fringes of the subluminal emission-see Fig. HOf a)). This 
method (Fig. HOf b). dotted curve) yields a value 2k = 0.90, i.e., an intensity proportional to 1/R 090 . 

For smaller values of R, a more sophisticated fitting method (Appendix; Eas. ll8l[TT)|l takes account of the beaming 
by using a Gaussian beam of angular width 2w, and also allowing for the presence of residual spherically-decaying 
waves close to the array. This function is shown as a solid curve in Fig. HOf b); the fit parameters are 2k = 1.05, 
2w = 4.18° and K — 16.6 (see Ea. I19|l. This implies that the non-spherically-decaying part of the radiation has an 
intensity that decays as 1/i? 105 and that it is confined to an angular width ~ 4°. The fact that K, which depends 
on the geometry and perfection of the source, is substantially larger than 1 implies that the dominant emission from 
the experimental machine is spherically-decaying, with a small component of non-spherically-decaying emission. 

The narrow angular width of the non-spherically-decaying part of the radiation is emphasised by Figs. HOf c) and 
(d), which show the detected power versus distance for <^h = and </>h = —10° (rather than the <j)n = —5° used in 
Fig HOf b)'). The power falls off more quickly at large distances in (c) and (d) than in (b), and this is reflected in 
the fits. For example, in (c), the simple model (dotted curve; Ea. I17|l yields 2k = 1.41 (with N constrained to be 
1.074 to match the fit in Fig. [THT a)). whereas the more sophisticated Eas. Hlsl and IT§I fsolid curve) give 2k — 1.47, 
2w = 6.69° and K = 17.3. Both fits suggest an intensity that is approximately oc 1/R 5 and the latter is suggestive 
of less beaming of the radiation. Similar considerations apply to the fits to the data in Fig. ITUT dl . 



C. Determination of the rate at which power falls off with distance using intensity ratios 



An alternative, model independent way of determining the exponent 2k is to divide the power detected along the 
cusp by that detected when the source is run subluminally under similar conditions. Such a comparison can only 
be made beyond the region containing the interference fringes (i.e. R > 150 m), because (as mentioned above) the 
cusp has a distinctly "beamed" character, whereas the radiation from the subluminal source is much more isotropic. 
This comparison is shown in Figs. lUf a) and (b) for data taken on two separate occasions with differing weather 
conditions. Each data set represents the detected power p cus P along the expected cusp direction with the source 
running super luminally (v/c = 1.064), divided by the detected power P sub with the source running subluminally 
(v/c = 0.875). In both subluminal and superluminal runs, all other experimental parameters were identical. 

The data were fitted to a function 

Tjcusp 

CR" (6) 



^sub 

where fi and C are fit parameters. The fits yielded fj, » 1, i.e. a straight line through the origin, implying that 
the power ratio is roughly proportional to R (Figs. UlT a'). (b)). This result can be explained as follows. Once the 
interference fringes have been passed, 'p sub will obey an approximate relationship 

nsuh 

V**>( R ) « __ (7) 

where C sub is a constant or a very slowly- varying function of R. There are many precedents for the use of such power 
laws in the empirical functions [Fa. l2ll Y22\ for describing radio transmission (e.g. Eq. 0}. The straight lines through 
the origin seen in Figs. Illf aUb - ) will only occur if the power p cus P follows the relationship 

/"»cusp 

^ CUSP -^, (8) 

so that 

'ir>cusp /~>cusp 

T}sub /^sub v ' 

This again suggests that the power in the non-spherically-decaying component of the radiation would decay in free 
space more slowly with distance R than 1/R 2 . 

Once again, the narrowness of the region occupied by the cusps can be demonstrated by shifting 0h by only 5° 
( Figs. nTT c) -(d)). The power ratio increases more slowly with distance, and the data can no longer be fitted by a 
straight line through the origin; instead fj, = 0.62 in (c) and 0.52 in (d) (see Eq. |HJl. 



D. Summary 



The various methods of characterizing the rate at which the detected power falls off with R shown in Figs. 1101 and 
llll suggest that the intensity of the non-spherically-decaying component is roughly proportional to 1/R, in agreement 
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with theoretical predictions |a,ll3- Moreover, the relatively small angular range over which this 1/R dependence is 
observed (compare Figs. HUf b) and (c)) emphasises the narrowness of the region occupied by the cusps. 

Note that in the present experiments, the Fresnel distance (Rayleigh range) is of the order of a few meters. The 
maximum length of the array that can be projected onto a plane perpendicular to the propagation direction is 
I = 1.74 m and the shortest wavelength used is A = 0.50 m. These together yield a maximum Fresnel distance 
i?F = l 2 /A = 6.0 m. There is no parameter other than the Fresnel distance to limit the range of validity of the above 
results, and these, as we have seen, have been tested for R ^> Rp. 

V. CONCLUSION 

In summary, we have constructed and tested an experimental implementation of an oscillating, superluminal po- 
larization current distribution that undergoes centripetal acceleration. Theoretical treatments |l2l Il8j lead one to 
expect that the radiation emitted from each volume element of such a polarization current will comprise a Cerenkov- 
like envelope with two sheets that meet along a cusp. The emission from the experimental machine has been found to 
be in good agreement with these expectations, the combined effect of the volume elements leading to tightly-defined 
beams of a well-defined geometry, determined by the source speed and trajectory. In addition, over a restricted range 
of angles, we detect the presence of cusps in the emitted radiation. These are due to the detection over a short time 
period (in the laboratory frame) of radiation emitted over a considerably longer period of source time. Consequently, 
the intensity of the radiation at these angles was observed to decline more slowly with increasing distance from the 
source than would the emission from a conventional antenna. The angular distribution of the emitted radiation and 
the properties associated with the cusp are in good quantitative agreement with theoretical models of superluminal 
sources once the effect of reflections from the earth's surface are taken into account. In particular, the prediction that 
the beaming and the slow decay should extend into the far-zone has been tested to several hundred Fresnel distances. 

The above results have implications not only for the interpretation of observational data on pulsars, which originally 
motivated the theoretical investigation of the present emission process |y, [Hj , but also for a diverse set of disciplines 
and technologies in which these findings have potential practical applications [26j . One such technology is long-range 
and high-bandwidth telecommunications. 

Since all of the processes involved in generating the polarization current distribution are linear, and the direction 
of the emitted "beam" only depends on the source speed, it should be simple to use an array to broadcast several 
beams in different directions simultaneously. This could be achieved by applying voltages which are sums of terms 
like Eq. [31 each with a different speed, to the electrodes of an array. The beaming is frequency independent, and so an 
array could potentially be used to transmit multiple directed beams at several frequencies. (The very small side lobes, 
especially in free space, are a further characteristic of the beaming produced by the present emission mechanism that 
facilitates this type of transmission.) 

Of particular practical significance are the properties of the non-spherically-decaying component of the radiation 
that we have observed: the fact that its intensity decays like 1/R, rather than 1/R 2 , with the distance R from its 
source, and that its beam-width is less than ~ 5° at ~ 500 m. This radiation mechanism could thus be used, in 
principle, to convey electronic data over very long distances with very much lower attenuation than is possible with 
existing sources. In the case of ground-to-satellite communications, for example, the above properties imply that 
either far fewer satellites would be required for the same bandwidth or each satellite could handle a much wider range 
of signals for the same power output. 

The excellent correspondence between the experimental data and theoretical predictions supports the proposal that 
polarization current distributions act as superluminal sources of radiation, and that they can be used for studying 
and utilizing the novel effects associated with superluminal electrodynamics. The apparatus described in this paper 
represents just one possible geometry and excitation scheme. In principle, a whole class of superluminal sources are 
possible, corresponding, for example, to a rectilinearly accelerated polarization current @ (rather than a centripetally 
accelerated one), or to a medium that is polarized by impinging laser beams (rather than by applied voltages), and 
so on. It is hoped that the publication of this paper will stimulate further work along these lines. 
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VII. APPENDIX: PROCEDURES USED IN MODELING AND FITTING THE EXPERIMENTAL DATA 

Introduction 

This Appendix contains the procedures used to model and fit the experimental data. It is necessary to be rather 
rigorous in treating the radiation from a superluminal source and the effects of interference from the ground, and so 
we have described the mathematical techniques in some detail, a level of detail that would impede the discussion of 
the experimental data if it were included in the main body of the paper. Moreover, some effects described in the 
Appendix, such as interference caused by reflections from the ground, feature in the interpretation in all sections of 
the paper. It is therefore logical to collect all of the discussions of these phenomena in one place. 

The Appendix is divided into two main sections, treatment of interference due to reflections from the ground (see 
Fig. 0}, and use of the theory contained in Refs. [12, [13 to model the emissions of the experimental array (and in 
particular the beaming of the radiation) . Finally a shorter section deals with the robustness of the (small number of) 
fitting parameters used to model the data. 

Interference caused by reflections from the ground 

It is well known from the telecommunications literature [l4l H7H21I that interference is a significant complication 
in the propagation of radio waves close to the ground; the detected signal is the result of the interference between 
the direct and the reflected waves, i.e. between the waves that propagate directly from the source to the detector 
and the waves that reach the director after reflection from the ground. In the current experiments, this effect has 
been characterized by observing the emission from both a small dipole aerial and the array of the source driven at a 
subluminal velocity. 

Simple treatment of the effects of interference on the radiation of a dipole aerial 

We treat the dipole aerial as a point source at x = 0, y — 0, z — h A emitting radiation of wavenumber k in the 
presence of a medium, index of refraction A, occupying z < 0. The detector is assumed to be at the height hp and a 
distance R away from the source (see Fig. 4(a)). 

Had it not been deflected by the ground, the reflected ray that passes through the observation point P would have 
passed through the mirror image P' of P with respect to the ground. Thus the angle a that the reflected ray makes 
with the normal to the ground e g , and the distance i? r that the reflected wave traverses to reach the detector are 
respectively given by 

a = arctan{[i? 2 - {h B - h A ) 2 ]?/(h D + h A )}, (10) 

and 

Rr = {R 2 +4h A h B )i (11) 

(see Fig. 4(a)). When R is large enough for the curvature of the wave fronts to be negligible, the reflected and the 
incident waves are related via the Fresnel coefficients: 

Ereflectcd = ?_L Evident! (12) 

or 



-^reflected <7|| Bi nc iclent 



(13) 
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(see ) . The vectors E and B in these relations are the electric and magnetic fields of the incident and reflected 
waves and 



(TV 2 



sin a 2 



cos a + (N 2 — sin a) 2 



(14) 



N 2 cos a - (TV 2 



sin a ) 2 



TV 2 cos a + (TV 2 — sin 2 a) 2 



(15) 



are the Fresnel coefficients for Evident perpendicular or Evident parallel to the plane of incidence, respectively; 
the magnetic permeability of the ground is here set equal to unity. (Note that the index of refraction N and so 
the coefficients q± and qu are in general complex. However, for waves of frequency /, the imaginary parts of these 
quantities are negligible when 2a/ (e.f) <C 1, where a and e are the conductivity and the dielectric constant of the 
reflecting mt 
this paper.) 

Hence, in the case of the spherically-decaying radiation that is generated by a point-like aerial in the vicinity of the 
ground, the observed power at the detector is given, to within a constant of proportionality, by 



reflecting medium, respectively 



The latter condition appears to encompass all of the experiments described in 



exp(i/ci?) 



exp(i/ci? r ) 



(16) 



where the subscripts in Vf, and V\ 



R i? r 

refer to the common polarizations of the incident and reflected waves: parallel 
or perpendicular to the plane of incidence. By measuring Vf< or V\ as a function of the distance R and fitting the 
resulting data to the above equation, it is possible to determine N for the ground (i.e. the airfield runway) under a 
particular set of weather conditions (see Fig. 4(b)). This is then used in Eqs. 1171 and ITD1 (see below) to model the 
emission of the superluminal source. 

The effects of interference on the non-spherically-decaying component of the radiation from a super- 
luminal source 

The radiation by a superluminal source consists of two components: a spherically-decaying component whose intensity 
has the dependence R~ 2 on R (as does any other radiation) and a nonsphcrically-decaying component. To treat the 
latter component in the presence of reflection from the ground, we rewrite Eg. 1161 as follows: 



exp(ifci?) 



R K 



■<?iu - 



exp(ifci? r ) 



R T K 



(17) 



In other words, the exponent n, describing the rate at which the amplitude of each component (direct and reflected) 
diminishes with increasing R, has become an adjustable parameter. The theoretically predicted value of k for a 
centripetally accelerated superluminal source is h (see 0,01)- 

Ea. 1171 will only be valid if the amplitudes of the radiation traveling via the direct and reflected paths are similar. 
However, in the present experiment, the beam-width of the nonspherically-decaying component is appreciably (by a 
factor of 10 or more) smaller than that of the spherically-decaying component. This results in a substantial difference 
between the amplitudes of the direct and reflected rays when the detector is close to the source. (It is this phenomenon 
that is responsible for suppressing the interference fringes in Figs. HOf bh (c) and (d).) Unless R is large, therefore, 
the Fresnel coefficients in the above expression should be modified to take account of beaming. When modeling the 
data for ranges R > 150 m, we have consequently replaced the coefficients q\\ j_ by 



Q\\,± = Q\\,± exp[-(|7r - a) 2 /w 2 }, 



(18) 



where a is the angle between an incident ray and the normal to the reflecting surface (see Fig. 1121) and w represents 
the Gaussian half- width of the beam. The exponential factor multiplying the Fresnel coefficients greatly reduces the 
amplitudes of the reflected waves for which a differs from 7r/2 by more than w. 

As R tends to a value much larger than c/ui along this narrow beam, the nonspherically-decaying component of the 
radiation dominates over the conventional component. At shorter distances from the source, however, the detected 
power receives contributions from both of these components. If the fitting includes data for relatively small values of 
R (see e.g. the solid curves in Figs. llOf b'l and (c)), both spherically- and non-spherically-decaying contributions must 
be combined as follows: 



Pi 



l.-L 
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exp[i(fci? + vr/4)] exp(ifci?) 



R K 
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Qw 



exp[i(fci? r + 7r/4)] exp[i(fci? 



R? 



K- 
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(19) 
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FIG. 12: Relationship between the positions of the source, the detector and the notional "image" detector. 



Here, the constant K stands for the relative strengths of the sources of the two components, and 7r/4 is the phase 
difference betweem them (see [l2l Il3| ). 

Extension to a volume source: description of the direct and reflected rays in terms of the experimentally 
measured coordinates 

The modeling of the data is based on the same coordinates as those used in the theory papers [T2I IT^ | , shown in 
Fig. EHc,d). In these coordinates, the source consists of the arc — 7r/36 < p < 7r/36 of a circular strip, centered on 
the origin, with the radius a — 10.025 m and the thickness Ar = 0.05 m. The source velocity points in the direction 
of increasing ip (see Fig. EUd); its angular velocity points in the direction of positive z). Since the position of the 
source in FigjiJa) is such that its velocity points from the left to the right, the spherical coordinates (Rp,9p,(fp) 
of Refs. [I2j, Il3| have to be based, in the present case, on a Cartesian system whose z-axis points downward, as in 
Fig- Etc). (The midpoint of the source arc lies on the x-axis of this underlying Cartesian system; see Fig.^c).) 

For comparing the data with the theoretical predictions, we need to express (Rp, Op, ipp) in terms of the experi- 
mentally measured coordinates (R,-dv,4>v) and (R, #h,</>h) of Figs.|2Ia) and|3fb). 

The Cartesian system (x',y',z r ) corresponding to the spherical coordinates (R, "dy, cf>v) of Fig. EUa) differs from 
the Cartesian system (x, y, z) corresponding to the spherical coordinates (Rp, 9p, pp) of Fig. |3Jc) in that its origin is 
shifted by the array radius a along the positive x-axis, its x and y axes are interchanged, and its z-axis is inverted: 

x' = y, y = x — a, z' = —z. (20) 

The Cartesian and polar coordinates of the observation point P in the reference frame shown in Fig. I^a) are related 
by 

x' = R cos -#v cos 4>v ! y' = R cos $v sin cf>y , z' = Rsmd. (21) 
These, together with the corresponding coordinates of P in the reference frame shown in Figs.|2Ic), (d), 

x = Rp sin $p cos (pp , y — Rp sin 6p sin ipp , z = Rcos0p, (22) 
yield the first required transformation: 

R P = (R 2 +a 2 + 2aRcosd v smcf> v ) 1/2 , 



9p = arccos ( — — — sin^v) , 
V Rp J 



Pp = — 3n/2 — arctan ( tanc/iy H t- —J. (23) 

V it COS Vy COS 0y ' 

In Fig.[3fb), the position vector of the observation point P is projected onto the (x' , z') plane of the shifted Cartesian 
system (x',y', z'), so that 

x = R cos 0h cos i9h , y — R sin 0h , z = R cos 0h sin -dpi . (24) 
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Rewriting these expressions for (x',y',z') in terms of (Rp, Op, i^p), via [x, y, z), and solving for Rp, 8p and (pp, we 
obtain the second required transformation: 

R P = {R 2 + a 2 + 2ai?sin0 H ) 1/2 , 



( - -^-cos^h sini? H ^ 



tp P = -3?r/2 - arctan ( tan — + —). (25) 

V cos vh R cos 0h cos »h ' 

Having expressed the position vector xp = (Rp, Op, (pp) of the detector in terms of the experimentally measured 
coordinates, we now know how to determine the characteristics of the direct ray xp — x that links each source element 
x = (r, ip, z) of the array to the observation point P. The next step is to determine the corresponding characteristics of 
the reflected rays that reach P from various volume elements of the array. As we have already seen in the treatment 
of a simple dipole aerial, both the angle a that a given reflected ray makes with the normal to the ground, e gi and 
the distance R r that this wave traverses to reach the detector at P (see Fig. I12|) can be easily inferred from the 
position vector xp/ of the image P' of P with respect to the ground. We proceed to determine xp< in terms both of 
the coordinates (R, $v,4>v) and of the coordinates (R, $h, 4*h)- 

Suppose that the plane of the array is vertical, as in Fig.Eta), and that the upward-pointing normal to the ground 
at the point of reflection, e g , makes the angle <fio with the positive x-axis of the (x, y, z) frame (whose z-axis is in this 
case parallel to the ground), i.e. 

e g = cos 4>oe x — sm4>Qe y . (26) 

(For a smooth and perfectly horizontal ground, (f>o is related to 4>v via (f>o — 4>y = arcsin[(/iA — hp,)/R\.) Then the 
position vector, in the (x, y, z) frame, of the image P' of P across the earth's surface can be written as 

xp/ = R + ae x — 2hp,e g , (27) 

where R is the position vector of P in the (x' , y' , z') frame, the frame whose origin lies at the midpoint of the array 
(see Figs. EJc) and El). 

The following expressions for xp/ and R in terms of their components in the (x,y,z) and the (x',y',z') frames, 

xp/ = Rpi sm9pi(cosipp>e x + sinippie y ) + Rpi cos9p/e z , (28) 

R = i?cos^v(cos(/)vea;' + sin <fiye y i) + Rsmi)ye Z ' , (29) 

together with the relationships 

&x' ^yi ^y' &z' &z (^0) 

between the base vectors of these two frames (see Eq. 12011 therefore yield 

Rp< = [R 2 + a 2 + (2h D ) 2 - 4aft, D cos 4> + 2aR cos z? v sin V - 4ft, D i?cos # v sin(</> v - 4>o)} 1/2 , 



o m , /i?cosz? v sin(/>v + a - 2/i D cos0 o \ /Q1 , 
wp, — — 37T/2 — arctan — ; — : — ■ , (31) 

V R COS Vy COS <py + 2hp, sm (fio / 

for the spherical coordinates of the image point P' in terms of the angles (^V) <t>\r) of Fig. GJa). 

To determine xp/ in terms of the coordinates (R, i)h, of Fig.^b), suppose that the upward-pointing normal to 
the ground, at the point of reflection, makes the angle i?o with the negative z-axis of the (x, y, z) frame (whose x-axis 
is in this case parallel to the ground), i.e. 

e g = — sin^oGy — cosi^oe^. (32) 
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(Here, too, $q is related to i?h via i?o — = arcsin[(ft,A — hp>)]/R when the ground is smooth and horizontal.) 
Inserting this alternative expression for e g in Eg. 1271 and using Eqs. [2111201 an d the following version of Eg. 1241 

R = R cos <^>h (cos i?h6i' + sin $h£V ) + R sin 4>ne y > , (33) 

we obtain 

R p/ = [R 2 + a 2 + (2h D ) 2 + 2aRsin(j) ii - 4h D Rcos<f> H sin(tf H - i?o)] 1/2 , 



n ( R cos 0h sm $h — 2/id cos $o s 

^p/ = arccos ( 

V Roi , 



( R sin dm + a \ , 

w P / = -3tt/2 - arctan ( ■ „ — - ) , (34) 

Y ' Vi?cos</> H cos?? H + 2/i D sin?V 

for the spherical coordinates of the image point P', with respect to the (x, y, z) frame in terms of the angles (Ah, 4>h) 
ofFig.l^b). 



Theoretical treatment of the spherically-decaying component of the emission from the experimental array 

We now describe the use of the theory of Refs. 0, ^1 to model the emissions of the experimental array. Our 
purpose is to predict the beaming of the spherically-decaying component of the radiation; relevant data are shown in 
Figs. 15191 (The detailed behavior of the nonspherically-decaying component will be left for a future work.) The same 
mathematical notation as that in Refs. |l2l Il3j| is used; for consistency with the original papers, we employ Gaussian 
electromagnetic units. 

We commence with the following far-field expressions for the electromagnetic fields E and B in the absence of 
boundaries [Eqn. (14) of Ref. [l^]: 

E.i /d 3 ,dt ^ p -f- |Xp - X|/C) Ax(nxg) ; B = n x E. (35) 
c l J |xp — x| V at J 

Here tp is the time at which the field is measured at the observation point P, and x and t are the position of a 
source element and the retarded time, respectively; n = xp/|xp| is a unit vector along the direction linking the source 
that is localized in the vicinity of the origin and the observer at P. For the polarization P described in Eq. [21 whose 
distribution pattern both rotates (with an angular frequency uS) and oscillates (with a frequency f2), the source term 
in the above expression is related to the cylindrical components (s r , s v ,s z ) of the vector s that designates the direction 
of the current density j = dP / dt by 

nx(nx<9j/<9i) = \lo 2 X^ =(tl± y 2 exp[—i(p<p - VLip / uj)]{[s v cos9 P sin(<£ - tp P ) 
—s r cosOp cos(<p — tpp) + s z sin#p]ej_ + [s r sin(<^ — (fp) 

+s<p cos((p — <pp)]e\\ } + {m — > — m, — > — Q}, (36) 

where the cylindrical coordinates (r, ip, z) and (rp, tpp, zp) of x and xp are with respect to the (x, y, z) frame shown 
in Fig. Etc), tp = ip — tot, fi± = (f2/w) ± m, ey = e Vp (which is parallel to the plane of rotation) and ei = nx en are 
a pair of unit vectors normal to the radiation direction n, and the symbol {to — > —to, f2 — > —0} designates a term 
exactly like the one preceding it but in which m and £1 are everywhere replaced by —to and — fi, respectively. This is 
Eq. (18) of Ref. [13. 

We are interested in the context of the present measurements in the radiation that is polarized parallel to the plane 
of the array; note also that the component s v is zero in the experimental array described in the current paper. Hence, 
the relevant component of the source term in Eo. 1361 is 

e || • [n x (h x dj/dt)] = \s r si\\(p — </5p)Kj X)±(^ ^ tow) 2 exp{i[(f2 ± muj)t =p m<^]}|, (37) 

where denotes the real part of Z. The fact that the array has the form of an arc means that the domain of 

integration in Eq. [23 consists, in the present case, of the cylindrical strip —5° < tp < 5°, 10.00 m< r < 10.05 m, 
—0.005 m < z < 0.005 m, that bounds the volume occupied by the array (see Fig. Since the r and z dimensions 
of this volume are appreciably shorter than the wavelengths at which the radiation is generated and measured in 
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the present experiments (0.50 m and 0.59 m), the integration with respect to these two coordinates may be omitted 
when evaluating the integral in Eq. 1351 without introducing any significant error. (It must be emphasized that this 
approximation may only be used in the case of the spherically-decaying component of the emission; it is necessary 
to include the full volume of the source when calculating the nonspherically-decaying component. For the current 
purposes modeling the spherically-decaying data of Figs. l5T9T the approximation is adequate.) 

Inserting Eq.|23in Ea. 1351 and performing the trivial integration with respect to t, we find that the electric field of 
the radiation that is detected at either the frequency /+ = Q + ttllu = n + uj or the frequency /_ = \Q — j7ilu\ = \n—\ui 
(see Tabled Jljill3) is given by 



r r /36 - 

E(/±) oc M< exp(— m±tdtp) / dip sin(ip — ipp) exp[i(n±R± rrup^/Rye^p (38) 

L J-tt/36 J 



fir/36 
'-tt/36 

to within a constant of proportionality depending on the source strength. Here, 

R = |x P - x|w/c = [Rp + f 2 - 2rRpsm9 P cos(ip ~ fp)]K (39) 

in which f = ru/c is the source speed in units of c, Rp = Rpuj/c — (rp 2 + zp 2 )^ , and 9p = arccos(zp/i?p) [i.e. 
(-Rp, 9p,ipp) comprise the spherical coordinates of P in units of the light-cylinder radius c/ui]. 

The expression in Eq.[2Hl yields the electric-field vector of the wave that propagates from the source to the detector 
directly. The corresponding expression for the electric-field vector of the reflected wave is given by a similar expression 
in which the coordinates (Rp, 9p,pp) of P are replaced by the coordinates (R T , 9p> , ppi) of the image P' of P and the 
integral is multiplied by one of the Fresnel coefficients 511 or q± (see the preceding subsections): 



E r (/±) oc 3J< gii j_ exp(— in±ujtp) / dtp sin(<^ — <ppt) exp[i(n±-Rr ± m(p)]/R T >e T , (40) 

1 J-tt/36 ' ' 



where 

R T = [Rp, + f 2 - 2rRp, sinf5 P , cos(<^ - cp P /)]*. (41) 



The reflected radiation is polarized along 



(42) 



if E(/±) is normal to the plane of incidence, and along 

G rw — 2(eg • e^, p/ )eg ^^ p/ (^^) 

if E(/±) is parallel to the plane of incidence, where e g is the upward-pointing unit vector perpendicular to the ground 
(see Ecis.1121 and !13|) . The coordinates (Rp>, 0p>, <pp*) have the values derived in Eg. 1311 when the array is positioned 
as in Fig. EJa), and those derived in Eg. 1341 when the array is positioned as in Fig. 01b). The values of q» and q± 
appear in Egs.ll4landll5l 

Next, we need to find the power that is received at P by an observer whose detector only picks up radiation of a 
given polarization, e.g. a polarization along the unit vector p. This can be done by adding the components p • E(/±) 
and p-E r (/±) of the electric-field vectors of the direct and the reflected waves and averaging the square of the resulting 
expression over the observation time tp. The result is 



P\\,x(f±) « 



P ■ <W dip sin(<p - ip P ) exp[i(n ± ^ ± m<p)]/R 



■ e rw Jli% 6 dip sm(p - ip.p>) exp[i(n ± _R r ± mip)]/R r 



(44) 



a general result that applies to both frequencies, irrespective of whether the radiation is polarized parallel or perpen- 
dicular to the plane of incidence. 

When the array is positioned as in Fig. |3fb), the generated radiation is polarized horizontally, i.e. in a direction 
normal to the plane of incidence. In this case, p • e^p and p • e rw in the above equation are equal and so can be factored 
out and absorbed in the constant of proportionality. But when the plane of the array is vertical, as in Fig.[3Ja), the 
electric-field vectors of the generated and reflected waves, though both in the plane of incidence, are no longer parallel 
to each other. In such experiments, we have used a receiving antenna whose polarization is vertical, i.e. for which 
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p = e g . According to Eas. 1261 and 1431 therefore, the factors in question have the values p • e^p = — sin(</?p + 0o) and 
p • e rw = — sin(yjp/ + cj>o) in these experiments. 

Note that the emission is beamed when the phase of the rapidly oscillating exponential in the integrand of Eq. |2H 
has a stationary point within the domain of integration, i.e. when 



d 

— [n±R ± m(<p - tpp) 



(45) 



vanishes for a value of p in — 7r/36 < tp < 7r/36. This stationary point would be degenerate, and so the beaming 
would be tighter, if 



_0 2 

dp 



-^[n±R ± m(ip - ipp)] 



vanishes at the same time. The above two conditions are satisfied simultaneously when 

1 - m 2 /(n±R P ) 2 



sin Op - m 2 /(n±Rp) 2 



(46) 



(47) 



and so the phase in question is stationary at p = pp + 2n — arccos[m 2 /(n± 2 ffp)] (cf. 0,E1)- The parameters in 
our experiments are chosen such that the value of r implied by Eq. 1471 for an observer in the far field, lies within the 
radial thickness 10.00 m < r < 10.05 m of the array. 



Adjustable parameters: effects of the array orientation and surface corrugations on the modeling of the 

beaming data 

There are no truly adjustable parameters in the analysis of the beaming data (Figs. I5I9|I in this paper, especially 
when R > 300 m. The speed of the source is constrained by the setting of the experimental machine and the refractive 
index of the ground is determined by ancillary experiments with the dipole aerial (see Fig.0J). The measured heights 
ft, A and hp, and distance R are also used as input parameters. The remaining uncertainties come from the difficulties 
involved with modeling the reflection of curved wave fronts and small inaccuracies in the angular coordinates. We 
treat the latter first. 

In modeling the experimental data of Figs. 15191 variations in the parameters f?o (see Ea. 13211 and 4>o (see Ea. 12611 . 
the angles marking the orientation of the array, have a significant effect on the relative phases of the rays traveling 
on the direct and reflected paths (Fig. H*2l. Ideally, with flat ground (i.e. the runway) and precision measurements of 
the orientation of the array, these angles would be completely constrained by the measured experimental coordinates 
(f>y and i?h respectively, plus the known h^, hp, and R. However, there are two problems with this approach: 

1. The measurement of the angular coordinate desribing the tilting of the array is subject to relatively large 
fractional errors, ~ ±0.5°; 

2. The runway surface used for the measurement has distinct corrugations corresponding to fractions of a degree; 
over long distances, these can result in non- negligible shifts in the apparent position of the reflected ray. 

These uncertainties were encompassed by making do or </>o (as appropriate) somewhat adjustable; adjustments of less 
than ~ 1° from the expected value (i.e. of the order of the uncertainty in the experimental angles) were in general 
sufficient to obtain very close agreement between model and data. 

A second cause of uncertainty is the use of Fresnel coefficients in dealing with reflection from the ground; Fresnel 
coefficients are of course derived for plane wave fronts 0] . As is well known in the modeling of radar beams ^7| > 
wave fronts close to a source have significant curvature, leading to an apparent loss of the reflected radiation due to 
"phase smearing" . In the present case, the effect becomes noticeable for R < 300 m. 

A detailed treatment of curved wave fronts is very complex, and beyond the scope of this paper. However, a number 
of empirical approximations may be used Following these precedents, the reflection coefficient in simulations of 
the beaming involving distances R < 300 m was scaled by a factor qo, representing the loss of signal on reflection due 
to the curved wave fronts. 
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